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Abstract 

The present paper studies the escape mechanism in collinear three point mass sys- 
tems with small-range-repulsive/large- range-attractive pairwise-interaction. Specif- 
ically, we focus on systems with non-negative total energy. 

We show that on the zero energy level set, most of the orbits lead to binary escape 
configurations and the set of initial conditions leading to escape configurations where 
all three separations infinitely increase as t — > oo has zero Lebesque measure. We 
also give numerical evidence of the existence of a periodic orbit for the case when 
the two outer masses are equal. 

For positive energies, we prove that the set of initial conditions leading to escape 
configurations where all three separations infinitely increase as t — > oo has positive 
Lebesque measure. 
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1 Introduction 



In computational chemistry and molecular dynamics it is common to use 
mathematically simple models to describe the interaction between a pair of 
atoms or molecules (see, for instance, [4] and [7]). Set within the framework 
of the classical n-body problem, these models ignore quantum effects for the 
gain of significantly less computation time in numerical simulations. 

The most commonly encountered molecular potential is the Lennard- Jones 
potential (or the 6-12 potential), which in given by the sum of a repulsive term 
at present at short ranges due to overlapping electron orbitals, and a long range 
attractive term (the van der Waals force, or dispersion force). The history 
of the Lennard- Jones potential and the motivation of its use for modeling 
molecules can be found in [1]. The Lennard- Jones potential is often used to 
describe the properties of gases, it is particularly accurate for noble gas atoms 
and it is a good approximation for neutral atoms and molecules. 

In this paper we study a system formed by three mass points (atoms) con- 
fined to move on a line, with pairwise interaction given by an attractive- 
repusive potential of the form: 



W(r) : = - 




where e and r are given real positive parameters, and a, b G R such that 
4 < a < b — 1. Note that for a = 6 and b = 12 we retrieve the standard 
molecular Lennard- Jones potential. Thus one could think of our problem as a 
classical mechanics counterpart of a linear triatomic molecule. 

We analyze the dynamics arising at total energies h equal to and above 
the dissociation level h = 0. Specifically, we investigate escape or scattering 
solutions, with focus on 2 + 1 escape configurations (or binary formation), 
where two atoms remain close while the third separates, and 1 + 1 + 1 escape 
configurations, where all three inter-particle separations infinitely increase as 
t — > oo. 

Our approach relies on the choice of coordinates. Specifically, we use McGehee- 
type coordinates, where the size of the system, thought of as the radius of a 
sphere in the norm given by the moment of inertia, becomes one of the vari- 
ables. Further, using inversion diffeomorphic transformations and regularizing 
the equations of motion, we convert states of infinite size to states confined 
to so-called infinity manifold, that is, a compact manifold glued to the energy 
level set. In particular, the 2 + 1 and 1 + 1 + 1 escape configurations appear 
as invariant submanifolds (lines or equilibria) of the infinity manifold. While 
fictitious, the infinity manifold gives useful information about the real flow via 
the continuity of solutions with respect to initial data. 

Consequently, we analyze the flow on the infinity manifold in both h = 
and h > cases. We prove that on the zero energy level set, a zero Lebesque 
measure set of initial conditions leads to 1 + 1 + 1 escape configurations. For 
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mass-symmetric systems, where the outer atoms of the same kind, we demon- 
strate the existence a heteroclinic orbit connecting two 1 + 1 + 1 asymptotic 
states that is persistent for small mass perturbations. For strictly positive en- 
ergies we prove that the set of initial conditions leading to 1 + 1 + 1 escape 
configurations is of positive Lebesque measure. 

Our investigation is complemented by numerical explorations performed for 
the case h = 0. These suggest that for most of the initial conditions, the orbits 
lead to a 2 + 1 escape configurations. For mass-symmetric systems we signal 
the existence of a periodic orbit that is reminiscent of the Schubart orbit in 
the collinear three body problem of celestial mechanics (see [9], [10] and [11]). 
This is in a sense no surprise, since the Schubart orbit, a periodic solution 
with two collisions per period, becomes evident after the double collisions are 
regularized, taking the form of elastic bounces. We suspect that a "Schubart" 
periodic solution is present in any collinear three mass point system where the 
pairwise interaction is given by either an attractive potential regularizable at 
double collisions, or a small-range-repulsive/large-range-attractive potential. 
We also believe that the existence of such periodic orbits might be proven 
using variational methods as in [12]. We find these interesting questions that 
we intend to investigate elsewhere. 

The paper is organized as follows: in Section 2 we review briefly the equa- 
tions of motion and integrals for a collinear three mass point systems. We also 
classify the escape configurations and find an invariant subspace in the case 
of mass-symmetric systems. In the next section, we introduce McGehee coor- 
dinates and observe that for h > the flow does not admit any equilibrum 
states. In Section 4, escape dynamics is discussed for the case of zero total 
energy. We start by converting infinite size states into states at the origin. We 
further regularize double collisions, that, while physically impossible, appear 
as singularities of the vector field. The infinity manifold is defined, which, due 
to the regularizations of double collisions, shows up as a compact manifold. 
The flow on the resulting infinity manifold, its equilibria and invariant mani- 
folds are retrieved. This allows to determine the Lebesque measure of the set 
of initial conditions leading to 1 + 1 + 1 escape configurations. We conclude the 
section by presenting some numerically obtained results, including the periodic 
orbit mentioned above. In Sections 5, we focus on the case h > 0. After apply- 
ing an appropriate change of space and time coordinates, the infinity manifold 
Nh appears as a manifold with boundary, namely a cylinder with caps given 
by the (fictitious) double collisions. The equilibria of Nh are classified into six 
manifolds with boundary. The dimension of the stable and unstable manifolds 
of the equilibria corresponding to 1 + 1 + 1 escape configurations are then cal- 
culated using the concepts of overflowing and inflowing invariant manifolds. 
The main result is stated and proven in Theorem 5.5 together with Corollary 
5.6. 
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2 Equations of motion 



Consider a system formed by three material points (atoms) of masses mi, 
m 2 and m 3 , constrained to move on a line, with the pairwise interaction is 
given by a potential is of the form (1). Let the configuration of the system 
be described by q = (91,92,93) £ B^ 3 , where without loosing generality, we 
assume an ordering qi < q 2 < 93 at all times, and denote the momenta by 
P := (Pi>P2,P3) £ I^ 3 - The dynamics is then given by a canonical Hamiltonian 
system with the Hamiltonian: 

H(q,p) = ^p T M- 1 p + W(q), 

where M is the mass matrix defined by: 

M := [diag(mi)] . 
The potential W(q) consists of a two homogeneous terms sum: 

W(q) = -t'(q) + V'(q), 

f(q)- E V(q)== E jZ^S 

i<«<i<3 9? I i<«<i<3 |y» I 

with ocij and /3y real positive parameters and the exponents a and b such that 

4 < a < b — 1. In the case of the Lennard- Jones potential, when a = 6 and 

6 = 12, we have: 

^ := 4e ij (rg'') a > 0, Aj := ^(^')^ a > 0, 

where and are the characteristic parameter and the depth of the po- 
tential well, respectively, corresponding to the interaction of the atoms i and 
3- 

Since the potential has singularities at the collision points 
A := {q G 1R 3 : 91 = 92, 92 = 93, 9i = 93}, 
the associated Hamiltonian vector field 

q := M- l p (2) 

p ;= - W 

is well-defined on the set (R 3 — A) x ]R 3 . Standard results of the differential 
equations theory ensure, for given initial data (q(0),p(0)) G (1R 3 — A) x M 3 , 
the existence and uniqueness of an analytic solution (q(t),p(t)), defined on 
a maximal interval [0,t*), t* < 00. Analogously, one can work with intervals 
of the form (t*,0]. In case t* is finite, the solution is said to experience a 
singularity. 
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Since the system is conservative, along any solution (q(£),p(t)) we have: 



H(q(t), p(t)) = ^p T (t)M- 1 p(t) + W (q(t)) = const. =: h 



or, equivalent ly, 
1 



p T (t)M" 1 p(t) - 



1 


\qi(t) 


~<l2(t)\ a 




1 


\<h(t) 


-Q3(t)\ a 




1 


\qs(t) 


-<7i(*)l a 



«12 - 



"23 



"13 - 



012 


MO 


-q2(t)\ b - a 




023 


\Q2(t) 


-q 3 (t)\ b - a 




As 




-qi(t)\"- a 



h. 



(3) 



Remark 2.1 TTie ODE system (2) has no singular solutions. 

Indeed, assume that as t — > t* < oo ; potential W(q(£)) — > oo.. 5mce i/ie 
kinetic term is positive for all t, the left hand side of 3) approaches infinity, 
while the right hand side is constant. This is a contradiction. 

This is in agreement with physical intuition. The repulsive potential term is 
dominant at short range interactions and thus collisions, double or triple, are 
impossible. 

Using the linear momentum conservation, we reduce the dimension of the 
system. Thus, we restrict the configuration space to the linear subspace 

Q := {q G 1R 3 : m x q x + m 2 q 2 + m 3 q 3 = 0} 

and the momenta to 

P := {p Gi 3 : pi + P 2 + Ps = 0} 

and so the system (2) defines a vectorfield free of singularities on the four 
dimensional space (Q \ A) x P. 

Definition 2.2 A collinear atomic-like three mass point system where the 
outer mass points m\ and m 3 are of the same kind, i.e. mi = m 3 , au = a 23 
and P12 = P23 is called a mass-symmetric system. 

Remark 2.3 For mass-symmetric systems, the Hamiltonian becomes 



#(q,p) 



' (pi + pI) + t^tpI - 



1 



2m! 



2m 2 



m - q2\ a 



w - 5s 



"12 



P 



12 



I ft— a 



"12 - 
1 







12 



m - Qi\ 



\qi - q2 

"13 



ft— a 



(4) 



P: 



13 



192 - <?3 I 

In this case, a direct verification shows that the subspace 

CC := {q 2 =P2 = 0,qi + q3 = 0,Pi+P3 = 0} 



m - qi 



b—a 
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is invariant for (2). Physically, CC consists in motions having m 2 fixed in the 
centre of mass, while the outer mass points move symmetrically with respect 
to m 2 . 

We end this section by introducing two definitions. 

Definition 2.4 A solution (q(£),p(t)) of (2) is called a 2 + 1 escape config- 
uration (or binary formation), if either 

lim \qi(t) — 92 (t) | = oo and lim |g 2 (t) — q$(t)\ < oo, 

t— >oo t^co 



or 



lim \qi(t) — 92(01 < 00 and lim |g 3 (i) — 9 2 (i)| = 00. 



Definition 2.5 A solution (q(£), p(0) °/ caZZed a 1 + 1 + 1 escape con- 
figuration if 

lim \q 2 (t) - 9i(t)| = lim \q 3 (t) - q 2 (t)\ = 00. 

t— >oo t— »oo 



3 McGehee coordinates 

We introduce a new set of coordinates similar to those defined by McGehee 
in [81: 



r : = 



yV^q 



s = r x q 

v = r b/2 p T s 

u = r b/2 (p - t;Ms) 

and obtain the following system: 

f := rv (5) 
u = b - v 2 + u T M _1 u - ar b ~ a [/(s) + 6F(s) 
s = M _1 u 

u = ^ - f ^vu - (u T M _1 u)Ms + r fc ~ a (VC/(s) + aC/(s)Ms) - (W(s) + bV(s)Ms) 

Thus r 2 is the moment of inertial of the system and escape corresponds to 
r — > 00. Let 

S = {q G Q : r 2 = qMq = 1} (6) 

be the unit sphere in Q in the norm given by the moment of inertia. A point 
s G S defines a configuration for the system of particles and (r, s) G (0, 00) x S, 
(r, s) — > rs can be thought of as polar coordinates on Q\ A. Since u T s = 0, the 
velocity is broken into a rescaled radial component v and a rescaled tangential 
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component u, and the bundle 

T := {(q, p) G (Q \ A) x P | q G S, p T q = 0} 

can be thought of as the tangent bundle on S. Note that radial velocity v 
measures the rate of change of the size of the system, where as u defines the 
rate of change of the configuration. 
The map 

(0,oo)xMxT^(Q\A)xP 

(r,v,s,u) — > (q, p) := (rs, r~ 6 / 2 u + v Ms) 

is a real analytic diffeomorphism and in the new (r,v,s,u) coordinates the 
energy integral becomes: 

- (u T M- x u + v 2 ) - r b ~ a U{s) + V(s) = hr b . 

We further reduce the dimension of the system following ideas introduced 
by McGehee ([8]). For reader's convenience, we repeat here part of McGehee's 
construction. For more details, see [8]. We assume, without loosing any gen- 
erality, that the particles are ordered on the line such that q\ < q 2 < q 3 . The 
ordering is preserved along any orbit. Let 

5 = {s G S : si < s 2 < s 3 }, T = {(s, u) G T : s G S }, 

51 = {s G S : si < s 2 < s 3 }, Ti = {(s, u) G T : s G Si}. 

and let a = (ai,a 2 ,a 3 ) and b = (b 1 ,b 2 ,b 3 ) be the unique points on S with 
ax = a 2 < a 3 and b\ < b 2 = b 3 . Thus, the set Si is a closed interval with 
endpoints a and b, while So is the corresponding open interval. The endpoints 
a and b correspond to double collisions and are given by: 



ai = a 2 = 



m 3 



(m 1 + m 2 )(m l + m 2 + m 3 ) 



< a 3 



mi + m 2 



/ m 2 + m 3 

Oi = - J - ( , , r < h 

y mi (mi + m 2 + m 3 ) 



m 3 (m 1 + m 2 + m 3 ) ' 



mi 



\j (m 2 + m 3 )(mi + m 2 + m 3 ) 
Next we define a diffeomorphism between [—1,1] x R and Ti. Let 





l 


l 


l 







1 


-1 


M ■= 


l 


l 


l 


A 2 := 


-1 





1 




l 


l 


l 




1 


-1 






and 



A :-- 



1 



mi + m 2 + m 3 



mi m 2 m 3 
mi + m 2 + m 3 



1/2 



M~ l A 



2- 
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Then we have A T MA = M, Aq e Q for all q G Q, and 



q T MAq = 0, 
A 2 q = -q. 



(7) 
(8) 



Also, SL T A T Mh > 0. So, considering Q with the inner product induced by M, 
the matrix A is a rotation by 7r/2 in Q, a and b have unit length, and {a, Ah} 
is an orthonormal basis for Q. Note that a Mb is a constant depending only 
on the masses and that < a T Mb < 1. Let A be the smallest positive number 
such that 

cos 2A = a T M b 

and note that A is a constant depending on the masses only and < A < 7r/4. 
The map S : [— 1, 1] — > Si defined by: 



s -> s = S(s) : = (sin 2A)" 1 [(sin A(l - s))a + (sin A(l + s))b] 



(9) 



is a real analytic diffeomorphism and the following relation is true (see [8] for 
proof): 

S'(s) = XAS(s). (10) 
Let s := S~ 1 (s) and u := s T A T u. Then s £ [—1, 1] and u E R and we have a 
real analytic diffeomorphism 

[0, oojxlx [-1, 1] x R -> [0, oo) x R x Ti 

(r,t>, s,u) -> (r,t),S(s),MMAS(s)) . 

which, when restricted to [0, oo) x R x (— 1, 1) x R, is a diffeomorphism onto 
[0, oo) xlx T . Finally, the potential terms are transformed by the maps: 



£/:(-M) 
V :(-!,!) 



and are given by: 



U(s) = sin a 2A 



«12 



s -> C/( s ) : = f/(S(s)) 
s -> : = ^(S(s)). 



+ 



«23 



(&2-&1) ; j (i + s) a [(03-02)^1^- 

«13 



m»2 - oij (1+s) 



<• \ sin A 



sin A(l— s) 



;i- s )« 
(ii) 



V"(s) = sin fe 2A 



P 



12 



/7 7 \sinA(l+s) 
[t>2 - Ol) (1+s) 



+ 



23 



Pis 



([(62 - M^^] (1 + *) + [(«3 - a 2 )^^] (1 - s)) 



(12) 
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The potential terms U and V are strictly positive convex functions, with s = 
±1 boundary singularities of degree {—a) and (— 6), respectively. In physical 
terms, s = ±1 represent double collisions. For mass-symmetric systems, both 
U and V are symmetric under the transformation s — > (— s), with sq = a 
common critical point. 

In the new coordinates, system (5) becomes 

f := rv (13) 
v = + u 2 - ar b - a U(s) + bV(s) 

1 

s = —u 
A 

and the energy conservation reads: 

^ (V + t; 2 ) - r 6 - a t/(s) + V(s) = hr b . 

Since we considered a < b — 1, the above vectorfield is of differentiable class. 

Remark 3.1 There are no equilibria for positive energy levels. 

Indeed, equating the right hand side of system (13) to zero, and taking into 
account the conservation of energy, the following relations must be satisfied: 

u = v = 

(b-a)U(s) = (-h)br a 

r . b - a dU(s) _ dV(s) = o 
ds ds 

Since U(s) > 0, it follows that for h>0, the second equation above cannot be 
satisfied. 



4 Dynamics for zero energy 

4-1 Regularization of the motion 

The main focus of this section is on unbounded motions on the zero energy 
level. We begin by applying the change of coordinates p — 1/r together with 

x . = up (b-a)/2^ y .= vp (b-a)/2 
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Introducing dr/dt — p ( b a ^ 2 as a new time parametrization, system (13) 
becomes: 



p = -py (14) 
y = V + a; 2 - aU(s) + bp^V(s) 

1 

s = —x 
A 

x= Q-l) ^ + if/'( S )-i/-^'( S ), 

where we keep the same "dot" notation for the derivative d/dr, and where 
"prime" means derivation with respect to the s variable. The energy integral 
reads: 

1 -(x 2 + y 2 ) = U(s)-p b - a V(s). (15) 

The vectorfield (14) is well-defined for (p,y,s,x) G [0,oo) x I x (—1,1) x 
R. Escape orbits correspond now to orbits where lim p(r) = 0. In this case, 

the asymptotic escape configuration is given by lim s(r), and the asymptotic 

velocity of the system's size by lim y(r). 

From the energy relation (15), we deduce that p is bounded by the domain: 
V := | (p, y, s, x) G (0, oo) x R x (-1, 1) x R : p 6 " a < 



Since both U(s) and V"(s) are strictly positive, V is non-empty for any choice 
of the parameters. 

We now extend the vectorfield to include the double-collision endpoints s = 
±1 by using a rescaling of the radial coordinates p together with a Sundman 
type regularization. More precisely, we use two appropriate positive functions 
<f>(s) and /3(s) such that by applying the change of variables (p,y,s,x) — > 
(R, y, s, w) where 

p:=(/){s)R, x:=f3{s)w (16) 
together with an appropriate time reparametrization, 

- = — (17) 

da /5 2 (s)' 1 ' 

the vectorfield becomes well-defined for all (R, y, s, w) E [0, oo) xRx [—1,1] x 
R. An optimal choice for and (3 is given by: 

0( S ) = W, fts)= (1 _ 1 g2)a/2 , (18) 
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thus obtaining: 



V = (f - l) (1 - s 2 ) a {y 2 - 2U{s)) + bR b ~ a (l - s 2 ) b V(s) + 



R = 


-R 




fa 


y = 


u 






s = 






/ a 


w = 






u 



(l-s 2 ) a y- -2s(l -s^^w 
A 



(19) 



w = Q - (1 - s 2 )yw - ^as(l - s 



2\ S_i o 



+ i(l - s 2 )^ E/'(s) - \R b - a (l - s 2 ) 6+ t y'( s ). 
A A 

with the energy integral: 

l -w 2 + - s 2 ) a (y 2 - 2U(s)) + R b ~ a (l - s 2 ) b V(s) = 0. 

Thus, at the double-collision endpoints s = ±1 now we have 

,.„- , . /2sin2A\ a ay> 
lim 1 - s 2 ) a U(s) = ^— =: K_ < oo 

/ 2sin2A\ a «23 



(20) 



and 



lim (1 - s 2 ) a U(s) 

s— *1~ 



lim (l-s 2 ) b V(s) 
lim (1 - s 2 ) b V(s) = 

s— »1~ 



A / (a 3 -a 2 ) a 



/ 2sin2A^ 
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a ; (&2-&1) 6 



/ 2sin2A s 



'23 



A J (a 3 -a 2 ) b 



K + < 00 
L_ < 00 
L + < 00 



lim — (l-s 2 ) a L>(s) 
lim [ 4- (1 - s 2 ) fe V(s) 



s->l" yds 

hm (^a-* 2 ) 6 ^ 



+ 



It follows that system (19) it is at least differentiable for all (R,y,s,w) G 
[0, 00) xlx [—1,1] xl provided 4 < a < b — 1. The domain of motion for R 



is 



£>= <(R,y,s,w) G [0, 00) xlx [-1,1] xl : R b ~ a < G — s 2 ) a U{s) \ 

Since both (1 — s 2 ) a £/(s) and (1 — s 2 ) b V(s) are well-defined and continuous 
on [—1, 1], it follows that they attained their extrema. Thus R is finite at all 
times. 
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Remark 4.1 Under the above reparametrizations, escape definitions (2.4) 
and (2.5) read as follows: 

A 2 + 1 escape configuration, (or binary formation), is a solution of (19) such 
that lim R{cr) = and lim s(a) = ±1. 

A 1 + 1 + 1 escape configuration is a solution of (19) such that lim R(a) = 
and lim \s(a)\ < 1. 

(J — *oo 

Remark 4.2 For any escape solution, the asymptotic configuration and the 
asymptotic velocity of the system's size are given by lim s(a) and lim y(cr), 
respectively. 

4-2 Dynamics for zero energy 
From the energy integral we have: 



R b- a{l _ s yv(s) = (1 - S 2 ) a ( U( S ) - X / ) - ':>r~ > 



1 

2 



(21) 



which substituted into the second and the last equation of (19) leads to: 

1 



R = -R 



(1 - s 2 ) a y- -2s(l - s 2 )i~ l w 
A 

2"l) w 2 + (b-a-2)(l-s 2 r (U(s)-\y 2 



(22) 



w 



s — — (1 — s 2 ) 2 w 



A 



= (f-i) (i-*V+^(i - « 2 )S- ] 



+ A (1_a)T 



2as 



w 



The dynamics decouples since the last three equations may be solved inde- 
pendently (if the functions (y(-), s(-),w(-)) are determined, then R(-) is given 
by relation (21)). We thus focus on the reduced system: 



V = - (\ - l) ^ + (b - a - 2)(1 - s 2 T (U(s) - \y 2 



s = —(1 — s ) 2 w 
A 



(23) 



a _\ „ 2 1 



2\S-1 



«' = I 2 — -*-J — s )yw + 7^(1 — ) 5 



+ A (1 " S) 



:i ~ * 2 )H*) 
y ' V(s) \ 2 y 



2as 



w 
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where (y,s,w) Glx [—1,1] xl such that 

(l - s*T (u(s) - \y 2 ) - \w 2 > 0. 

Remark 4.3 The flow of the system (23) is symmetric under (y,s,w,a) 
(-y,s,-w,-a). 



The equilibria are found at the points P ± := (±y2U (s e ), s e , 0), where s e is 
the critical point of U (s). From relation (21), we deduce that at the equilibrium 
R = and thus P ± correspond to escape orbits. Also, since \s e \ ^ 1, the 
equilibria P ± correspond to 1 + 1 + 1 escape configurations. In particular, we 
have: 

Remark 4.4 For zero energy, in the case of 1 + 1 + 1 escape configurations, 
the asymptotic velocity of the system 's size is finite. 

The local behaviour near equilibria is determined using the linearization of 
(23) at P ± . A direct calculation shows that at P + the eigenvalues are given 
by 

£ 1 = -(l-sl) a y/2U(s e )<0 
and £ 2 < and £ 3 > 0, the roots of 



= 0. (24) 

Note that the eigenvector corresponding to £i is aligned to the y axis, 
where as the eigenvectors v& and v^ 3 are directions in the s — w plane tangent 
to M at P + . 

Since all eigenvalues have non-zero real part, the rest point at P + is hyper- 
bolic. We deduce that this equilibrium has a two dimensional stable manifold 
and a one dimensional unstable manifold. An analogous calculation shows that 
the stable manifold of the equilibrium at P~ is one dimensional, whereas the 
unstable manifold is two dimensional. So we have proved: 

Proposition 4.5 On the zero energy level, the flow of the collinear attractive- 
repulsive 3-particle problem, as defined by system (23), has two hyperbolic rest 
points, denoted P ± . Their stable manifold W s and unstable manifold W u have 
dimensions as follows: 

dim W U {P + ) = 1 dim W S (P + ) = 2 
dimW^(p-) = 2 dimW^(P^) = 1. 

Since the equilibria at P ± correspond to 1 + 1 + 1 dissociations, we have 
the following: 
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i \{i-sirpu{s e )^ 



- siy u d*U 
T 2 ~d& 



Corollary 4.6 In the linear attractive-repulsive 3-particle problem given by 
(2) with zero energy, the set of initial condition leading to 1 + 1 + 1 escape 
configurations has zero Lebesgue measure. 



The dynamics of (23) admits the manifold 




0}, (25) 



invariant. We call M the infinity manifold] it is a two dimensional bound- 
ary submanifold glued to the three-dimensional zero energy surface and it is 
formed by orbits with R = for all times. Thus, orbits on M represent ficti- 
tious motions where size of the system is infinite. The real flow is located in 
the region bounded by M containing the origin. By using the good differen- 
tiable properties at the boundaries and the continuity with respect to initial 
data, we can extract information about real orbits close to those on M. 
The dynamics on M is given by: 



It is immediate that the flow is decreasing along the y coordinate, since y 
is always negative. The boundaries {s = ±1} appear as invariant manifolds. 
They are four one dimensional boundary submanifolds glued to M, given by 



Since {s ± 1} correspond to fictitious double collision configurations, we call 
DC^ 1 the double collision manifolds. The dynamics on DC ± is given by, re- 
spectively: 



Orbits tending to DC^ 1 correspond to 2 + 1 escape configuration. It follows 
that in this case, y, the asymptotic velocity of the system's size, is unbounded. 
In other words, for a binary formation, the system is uniformly increasing in 
size. This is in contrast to the case of 1 + 1 + 1 escape configurations, where 
the asymptotic velocity of the system's size is finite. 

In Figure 1 the qualitative features of the flow on M is depicted for mass- 
symmetric systems. For all other systems, the picture is qualitatively the same, 



w 




DC :— {(y, s,w) elx [-1, l]xR : s = ±1, 



l -w 2 - K ± = 0} (26) 




w = 0. 
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Fig. 1. The zero-energy infinity manifold together with its flow and the heteroclinic 
orbit joining the 1 + 1 + 1 escape configurations P~ and P + . 

only the figure's "neck" is shifted to the left or right of s = 0, depending on 
where the critical point s e of U (s) is situated. Using the continuity of solutions 
with respect to initial data we deduce: 

Remark 4.7 There exists orbits connecting two different 2 + 1 escape config- 
urations. 

Remark 4.8 There exists orbits connecting 2 + 1 to 1 + 1 + 1 escape config- 
urations. 

For mass-symmetric systems, due to the symmetry of U and V, the flow of 
(23) is symmetric under the transformations 

(y,s,w) -»• (y, -s, -w). (27) 

In this case, there exits a heteroclinic orbit connecting the equilibria at P~ 
and P + . This orbit that may be determined by integrating 

V=(b-a) (U(0) - \y 2 ) . 

Proposition 4.9 For mass-symmetric systems, the unstable manifold W U (P~) 
and the stable manifold W S (P + ) intersect transversally along the heteroclinic 
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orbit that connects P~ to P + . 

Proof. Consider an orbit 71 := {(y(a), s(a),w(a)) : a E R} G W U (P~). 
Using the symmetries given by Remark 4.3 and relation (27), we have 

(- v /2C/(0),0,0)= lim (y(a),s(a),w(a)) 

v (J — > — O0 

= lim (y(a),-s(a),-w(a)) 

o^ — 00 

= hm (-y(-a),-s(-a),w(-a)) 

(J — ► — OO 

= lim (-y(<r),-s(a),w(a)) 
= - Hm (y(<r),s(<T),-w(<r)) 

and so 

lim (y(<7), 5(a), -«;(<t)) = (+^2U(0), 0, 0). 

Thus, 72 := {(y(o-), s(<j), -w{a)) : a ER} E W S {P+). 

Assume W U (P~) = W S (P + ). In particular, it follows that 71 G W S (P + ). 
Further, near P + , both 7 X and 72 must be tangent to the stable eigenspace of 
P + . Thus, near P + , r y 1 and 7 2 may be expressed as a linear combination of 
the stable eigenvectors and v% 2 . 

Now, Vfr is aligned to the y coordinate. Thus, close to P + , the s — w com- 
ponents of 71 and 72, that is (s(a),w(a)) and (s(cr), —w(cr)) , either are both 
zero, or are parallel to v^ 2 . In the first case we retrieve the heteroclinic orbit 
that connects P~ to P + . In the second case, we must have w(a) = 0, and, 
since no restriction was done over 71, this leads to the fact that the stable 
eigenvector v^ 2 must be parallel to the s axis. The last statement can be easily 
verified that is not true. So we conclude that W U (P~) and W S (P + ) intersect 
transversally along the heteroclinic orbit that connects P~ to P + . 

Corollary 4.10 The heteroclinic orbit P~ P + is stable under parameter per- 
turbation, and in particular, it will be present in the mass-asymmetric case for 
small mass perturbations. 

4-3 Numerical investigations 

We complete our study by presenting some preliminary numerically explo- 
rations. These show that most of the initial conditions lead to 2 + 1 escape 
configurations. A typical orbit emerges from an asymptotic 2 + 1 configura- 
tion, approaches a 1 + 1 + 1 configuration, and falls back into one of the two 
possible 2 + 1 escape configurations. 

We also noticed that in the mass-symmetric case, a periodic orbit (PO) 
appears surrounding the heteroclinic orbit P~P + . The PO is surrounded by 
its normal modes, orbits that generically fill in densely a cylindrical band 
around the PO - see Figure 2 - . The PO seems stable under perturbations of 
the parameters such as the ratio of the masses and the Lennard- Jones a and 
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-0.28 



-0.08 . 12 
s(t) 



Fig. 2. A periodic orbit (PO) for mass-symmetric molecules together with a 
surrounding orbit. The parameters are mi = 1 and rri2 = 2, ctyi = 1, 
fi\2 = 0.5, and = 2, /3i3 = 1, a = 6, b = 12, with initial conditions 
y(0) = 0,s(0) = 0,w(0) = 6.67 for the PO, represented by the black curve, and 
y(0) = 0, s(0) = 0,^(0) = 6.5 for the surrounding orbit. The vertical red segment 
depicts the heteroclinic orbit P~P + . 




Fig. 3. A periodic orbit (PO) for mass-symmetric molecules together with an escape 
orbit. The parameters and the initial conditions for the PO were taken as in the Fig- 
ure 2. The initial conditions for the escape orbit were y(0) = 0, s(0) = 0, w(0) = 6.9. 
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(3 parameters. 

As mentioned in the Introduction, the PO is reminiscent of the Schubart 
orbit discovered in the context of the collinear Newtonian three body prob- 
lem with regularized double collisions. We suspect that the Schubart orbit is 
present in any collinear mass-symmetric three body problem where the pair- 
wise interaction is given by either a potential regularizable at double collisions, 
or a small-range-repulsive/large-range-attractive potential. We intend to fully 
investigate this topic elsewhere. 

5 Dynamics for positive energy 

For h > we consider a different change of variables and rescaling of time 
p = 1/r, v = p b / 2 v, s = s, u = p b/2 u, dr = p b/2 dr]. 
System (13) now becomes 



v = u T M _1 u - ap a U(s) + bp b V(s) 

u' = -uv - u T M- 1 uMs + p a (VU(s) + aU(s)Ms) - p 6 (W(s) + aV(s)Ms). 
together with the energy relation 



As in Section 3, we reduce the dimension of the system by introducing the 
new variables 



We also use the potentials U and V introduced in equation (11) and (12). In 
these reduced coordinates system (28) reads 




pv 



(28) 




§ = S- 1 (s), u = s T A T u. 



—pv 



v' 



u 2 - ap a U{s) + bp b V(s) 



1 



(29) 




1 dU 1 „ h dV 



and the energy integral takes the form 
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We further regularize double collisions using a Sundman type regularization. 
Let 

p = (i-s 2 )R = e( s )R, 3^ = e(s) 

da 

Then the energy relation becomes 



(u 2 + v 2 ) - R a Q(s) a U(s) + R b Q(s) b V(s) = h. 



(30) 



Using the energy relation we can write the equations of motion as 



ft = - \ < ^-uR-QRv 
A as 

v' =6(5) 



u - aR a e(s) a U(s) + bR b Q(S) b V(s) 



1 



s' =-Q(s)u 
A 

u' — — Q(s)uv + 



A 



#»e(s) a+1 ^ -R b e(s) b+ld f 

as as 



(31) 



where ' now denotes differentiation with respect to a. These equations define 
a real analytic vectorfield on 91 := [0, oo) x R x [—1, 1] x R. For each h > let 



F fe = {(/?,{;, s,u) G 9t : -(fT + ^) - iTe(§) a f/(s) + J R f, 9(s) f, \/(5) = h}. (32) 

Z 



Since the gradient of expression (30) does not vanish on F^, F^ is a three 
dimensional real analytic submanifold of 91. F h can be written as a union of 



N h = {(R,v, s,u) em: R = 0,-(u 2 + v 2 ) = h} 

z 



(33) 



D h = {(R,v,s,u) e 9t : s = ±l,^(u 2 +v 2 )-ftQ(s) a U(s)+ftQ(s) b V(s) = h} 

(34) 



E h = {(p,v,s,u) E 91 : R > 0, -(M 2 + w 2 )- J R a 0(5) a t/(5) + J R b e(s) b V r (s) = /i}. 

z 

The sets JV^ and are invariant under the flow (31). We call Nh the infinity 
manifold and Dh the double collision manifold. The infinity manifold and 
the double collision manifold appear as a boundary submanifold glued to the 
energy surface Eh- Note that iVj, is homeomorphic to a cylinder in the three- 
dimensional space (v,s,u). 
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5. 1 The flow on the infinity manifold Nh 



The flow on the infinity manifold is given by: 

v' = Q(S)u 2 

§' = \e(s)u (35) 
A 

v! = —Q(s)uv. 

Note first that (35) the only restpoint are (v, s, u) = (±y/2h, s, 0), (v,s,u) = 
(v, 1, ±^/2h — v 2 ), and (v, s, u) = (v, —1, ±^/2h — v 2 ). Note that the equilibria 
(v,s,u) = (±y/2h,s,0) correspond to 1 + 1 + 1 escape configurations. We 
classify the equilibrium points in the following six manifolds with boundary 
as follows 

7V± = {(v,s,u) 6Rx [-1, 1] x R\(v,s,u) = (±V2h,s,0)} 
Df = {(v,s,u) 6Rx [-1,1] x R\(v,s, u) = (v, l,±V2h-v 2 )} 
Dt x = {{v,s,u) 6Rx [-1,1] x R\(v,s,u) = (v, -1, ±y/2h - v 2 )} 

Proposition 5.1 The flow on Nh is gradientlike with respect to the v coordi- 
nate. 

Proof. We show that v increases on every nonequilibrium solution. From the 
first of the equations in (35) it is clear that v' > 0, since s G [—1,1]. If v' — 
then either u = or Q(s) = 0. In both cases we are at an equilibrium point. 
Corollary 5.2 If we exclude the equilibrium point there are no closed or re- 
current orbits on Nh. 

Proof. If so, v would vanish identically on such orbits. 

As a consequence, all orbits of the flow must tend toward one of the equi- 
libria. 

To study the flow on Nh we introduce the new variable x £ [ — ^i 71 "] via 

u = v2/icos y . , 

(36) 

v = v2hsmx 

we find 

% = A. (37) 
ds 

The behaviour of the flow changes with A. See Figure 4 for a picture of the 
flow for a particular value of A 

5.2 The flow off the infinity manifold Nh 

We now want to study the flow off N h . In order to do that we will introduce 
the concept of overflowing and inflowing invariant manifolds and we will then 
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N£ (x = tt/2) 



Nh (x = -tt/2) 



s 

Fig. 4. The flow on the infinity manifold Nh for /i > 0. 

prove the existence of stable and unstable manifolds for the invariant manifolds 
Nf^ and we will find their dimension. 

Definition 5.3 Let (f> t be a smooth flow. M = MUdM is said to be overflow- 
ing invariant if for every p G M , <fit(p) G M for allt < and £/ie corresponding 
vector field is pointing strictly outward and is nonzero on dM. 

Let <f) t be a smooth flow. M = M U dM is said to be inflowing invariant if 
for every p G M , 4>t{p) £ M for allt>0 and the corresponding vector field is 
pointing strictly inward and is nonzero on dM. 

M = M U dM is said to be invariant if for every p G M , 4>t{p) £ ^ /o r a ^ 

f. 

Let <j) t be a smooth flow on a manifold C and suppose M is an overflowing 
invariant submanifold of C. Let P u C TC|^ be a subbundle that contains 
TM and is negatively invariant under the linearized flow D<p t - Let I C P u 
be any subbundle complementary to TM, and let J C TC be any subbundle 
complementary to P".Then the tangent bundle to C over M splits into three 
subbundles TM, I and J, i.e. TC|^ = TM © 7 © J (where © is the Whitney 
sum). 

Let II 7 , II J and I1 T be the projections onto I, J and TM respectively cor- 
responding to this splitting. We define generalized Lyapunov type numbers as 
follows: For any p G M, let 

A(p) = inf{6 G R" 
7(p) = inf{a G R' 

and, if j(p) < 1, define 

trip) = inf{s G R : D(f)_ t (p)Tl T \\ \\Il J D<fo(<f>_ t (p)) \\ s -> as t -> oo}. (39) 




|n J D0_ t (p)||6* -> as -> oo}, 
||n J D^(0_ t (p))||/a' -> as t -> oo}, 



(38) 
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More computable expressions for the Lyapunov type numbers can be derived 
from the expressions above, and have the following form: 



X(p) = lim WU 1 D^-tipW 11 
7 (p)= / lim||n J J D0 t (0_ t (p))|| 1 
. — \n\\D^ t (p)Il T \ 




*-oo-ln||n^D0 t (<Mp))|| 



where || • || is some matrix norm. X(p), 7<jo) and a(p) are independent of the 
choice of the metric. They also have some other nice properties, see [2] for 
more details. 

We have the following theorem: 
Theorem 5.4 Let <pt be a smooth flow on a manifold C and suppose M = MU 
dM is a manifold with boundary overflowing invariant under <p t , with P u C 
TC\m a subbundle containing M negatively invariant under the linearized flow 
D(p t . Then if \{p) < 1, 7(p) < 1 and a{p) < for all p e M there exists a 
smooth manifold W u overflowing invariant under 4>t such that W u contains 
M and is tangent to P u along M . 

Furthermore the manifold W u is permanent under small perturbations of 
the flow. 

For the proof of the Theorem we refer the reader to [2]. We remark that 
Theorem 5.4 can be applied to inflowing invariant manifolds. In that case, the 
generalized Lyapunov exponents are computed using the time reversed flow 
with the limits taken as t — > — oo, and the phrase "overflowing invariant" 
replaced with "inflowing invariant". Also P u and W u are replaced by P s and 
W s , with P s taken to be a positively invariant subbundle under the linearized 
flow D(f) t . 

We are primarily interested in applying Theorem 5.4 to the line segments 
of equilibria N^. We can now prove the following: 

Theorem 5.5 The manifold has a three-dimensional local stable mani- 
fold, and the manifold Nj^ has a three-dimensional unstable manifold. 
Proof. Let cf = (0, ±V2h, 5, 0) and denote by X : 9\ -> R 4 be the vector 
field associated to the system (31). We want to compute the eigenvalues of 
the Jacobian 



where x = (0, v,s, u) G F h . In particular we want to determine the Jacobian 
for x — cf. Note that DX(x) : 1R 4 — >■ M 4 leaves T x F h invariant, and J is the 



J '■ T x F h 



T x F h . 
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restriction of DX(x) to T x F h . Using system (31) we have 

'2su 



[l -~s 2 )v)a 



(2su 

DJi(x)(a,7,(T,x) = I -y 

DI 2 (s)(a, 7, <r, X) = (-2sm 2 )ct + (2m0(5)) 7 

1 25^ 
DX 3 (a:)(a,7,(T,x) = (jQ{s))x + (-- ^-u)<r 



(40) 



M 4 (a;)(a,7,<T,x) = (-Q(s)u)-f + (2suv)a + (-Q(s)v) X 
where x = (0,5, s,u). From the definition (32) of F h we see that 

T x F h = {VF h ■ (a, 7 , a, x ) = 0} = {(a, 7, a, x ) e R 4 : ^ + u X = 0}, 
and, in particular, T c ±F h = {(a, 7, a, x) £ : ±^^7 = 0}. 

s 

Thus, for (a, 7, a, X ) G TL-F/i, we have 



■7 («, 7, o", x) = DX(x) (a, 7, a, X ) 

'' (1 .s 2 )r)n 



+ 



A 

2s 



A «)a + (^e(s)) x 



ei + 
e 3 + 



u 2 



{-2su')a + (-2^e(s)) X 



e 2 



« 2 



(2suv)a + (©(s)y - 0(s)w)x 



e 4 



where ei,e2,e3,e4 are the elements of the standard basis of 1R 4 . A basis for 
T x Fh is given by the vectors ei, e 3 , e 4 Then the matrix for J in this basis is 









(2Siw)<7 (6(5)^-6(5)15) 





(ie(s)) 



In particular, J at x = cf is given by: 



-0(5)5 







i©(5)) 







- 0(5)6 



(41) 



Thus ex is an eigenvector with eigenvalue =F©(s)v2/i. Note that {e 3 ,e 4 } is a 
basis for T c ±N h . The characteristic equation for jf restricted to T c ±N h is 



-£(-£ T V2/i0(5)) = O. 

The above equation has roots £1 = 0, and £2,3 = tV2~KQ(s). The eigenvalues 
of J(cg) + are 0, with multiplicity one, and — ®{s)\/2h with multiplicity 2. 
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The eigenvalues of J(cg) are 0, with multiplicity one, and Q(s)y/2h with 
multiplicity 2. 

Let M ± = {(v,s,u) ERx [-a, a] x R\(v,s,u) = (±V2h,s,0)} with a < 1. 
M ± is a proper subset of iV^. We now discuss the hyperbolicity properties 
of M~ and compute the generalized Lyapunov type numbers. Similar compu- 
tations hold for M + . Let us denote with E s (cf) and E u (cj) the stable and 
unstable subspaces of the linearized equations, consider the following unions 

(J (E s (cf)), E" = (J (E*(cj)). 

a, a] a,a] 

Then we have 
If we define 

pu _ 

P s = E s 



then it should be clear that P s is a positively invariant subbundle and P u 
is a positively invariant subbundle. Since TM ± = M ± x (0,0, R, 0), we can 
write TF h \ M - = TM~ © 7 © J, where J = M~ x (0,0, R, 0) and I C P" (a 
subbundle complementary to TM~). I and J are the subbundles / and J in 
the geometrical set-up of Theorem 5.4. Hence the projections (in the basis 
ei, e 3 , e 4 ) are 










1 









1 


, n' = 





, n J = 












1 








Now we compute the generalized Lyapunov type numbers associated with 
P u in the context of Theorem 5.4. Using the expression of the Jacobian com- 
puted at x — cf (41) we have 



DM4) 



e -e(s)i,t q 






I»-e(!> 



e e{ ~ s)u 



Let || • || be the matrix norm defined by \\A\\ = (trace (AA 1 )) 1 ' 2 . Then 

||n'£><Mcr)|| = trace((n 7 £>^_ t ( C r))(n 7 £>^( C r))*) 1 /2 = ^-e^V**), 
and hence 

X(cj) = Rm(v / 2e- e ^^) 1 /* = e' @ ^^ < 1 
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for a < 1. A similar computation shows that, since <f>(c s ) = c s , ||n J D0 t (c s ) || = 
0, and thus 

7(cj)= Iim||n J D0 t (cr)|| 1 A = O <l. 

r — >oo 

Furthermore, since ||il J D0 t (cj)|| = ||i^0_ t (cj)n T ||=O, then cr(cj) can be eas- 
ily computed to be using equation (39). 

At this stage we have to apply Theorem 5.4. This cannot be done directly 
since the vector field is identically zero on dM~ and thus M _ is not overflowing 
invariant. This is a technical detail that can be dealt by using a C°° bump 
function to modify the vector field in a small neighborhood of dM~ as in 
Section 18 of [3], in Section 4.1 of [13] and in Proposition 2.3 of [14]. This 
concludes on the existence of a three-dimensional local unstable manifold W% c . 

Similarly one can apply a variation of Theorem 5.4 valid for inflowing in- 
variant manifolds to prove the existence of a three-dimensional local stable 
manifold Wf oc for M + . This concludes the proof. 

From the theorem above we have the following: 

Corollary 5.6 In the collinear attractive-repulsive 3-particle problem with 
positive energy given by (2), the set of initial condition leading to 1 — 1 — 1 
asymptotic escape configuration has positive Lebesgue measure. 
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